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Abstract 

In a new measurement on neutrino oscillation ^ I'e, the MiniBooNE Collaboration observes an 
excess of electron-like events at low energy and the phenomenon may demand an explanation which 
obviously is beyond the oscillation picuture. We propose that heavier neutrino 1^2 decaying into a 
lighter one vi via the transition process ^ z^g + ^ where X denotes any light products, could be 
a natural mechanism. The theoretical model we employ here is the unparticle scenario established 
by Georgi. We have studied two particular modes — > Ve+^ and z^e + i^e + '^e- Unfortunately, 
the number coming out from the computation is too small to explain the observation. Moreover, 
our results are consistent with the cosmology constraint on the neutrino lifetime and the theoretical 
estimation made by other groups, therefore we can conclude that even though neutrino decay seems 
plausible in this case, it indeed cannot be the source of the peak at lower energy observed by the 
MiniBooNE collaboration and there should be other mechanisms responsible for the phenomenon. 
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I. INTRODUCTION 



Recently, the MiniBooNE Collaboration reported its results of searching for z/^ — >■ z/g 
oscillations In the experiment, the energy spectrum has a peak centered at 700 MeV 
and extends to 3000 MeV. For the oscillation range 475 < < 1250 MeV where is 
the energy of the produced neutrino, no significant excess of events is found. This result 
excludes sizable appearance of z/g via two neutrino oscillation and disfavors the previous 
LSND measurement However, they observed that outside of the oscillation range, there 
is a clear peak of the electron-neutrino-like events (96±17±20 events) lying above background 
at 300 < E^ < 475 MeV. Although the origin of the excess is still under investigation, we may 
assume that they are indeed electron neutrinos at present. The beam is completely composed 
of z/^ and the oscillation can only produce with the same energy, therefore the observation 
would be a serious challenge to the present theories. Namely a reasonable explanation about 
the appearance of the low energy is needed. To answer this question, there are some 
interesting proposals, for example, in jl] the authors suggest a (3+2) neutrino oscillation 
scenario where two sterile neutrinos are introduced into the game to explain the MiniBooNE 
results and BodekQ] considered the internal bremsstrahlung as an alternative source of the 
excess z/g events. Instead, in this work, we are looking for possible mechanisms other than 
the neutrino oscillation, supposing that there are only standard model (SM) neutrinos. An 
explanation that z/^ may decay into + X where X denotes some possible hght products, 
seems reasonable. Definitely, neutrino decay must be realized via interactions beyond the 
SM. The possible candidates of X could be z/g + z/g, light bosons (for example axion etc. ) 
and the unparticle which we are going to explore in this work. 

In fact, the idea of neutrino decay is not new. It has been put forward by some authors 
[sl, [gI. The basic idea is to introduce a heavy, unstable neutrino (usually assuming a sterile 
one) which decays into light neutrino or antineutrino plus a scalar particle. The interactions 
between the scalar particle and neutrinos are described by a lepton flavor violating effective 
lagrangian which depends on the details of various new physics models. Instead, we suggest 
an alternative scenario, namely the heavier z/2 which is a mass eigenstate and a component 
of the flavor eigenstate z/^, is the constituents of the beam and decays into a light neutrino 
z/i and a scale- invariant unparticle proposed recently by Georgi [3]. 

It is well known that at very high energy scale, the unparticle physics contains the SM 
fields and a sector of Banks-Zaks field (defined in [sl) with a non-trivial infrared fix point. 
Below an energy scale An which is of order of TeV, the Banks-Zaks fields are matched onto 
a scale invariant unparticle sector. The unparticle is different from the ordinary particles as 
it has no mass since the mass term breaks the scale invariance, but the Lorentz-invariant 
four-momentum square needs not to be zero, > 0. The scale dimension of unparticle is 
in general fractional rather than an integral number (the dimension for a fermion is half- 
integers). This special characteristic brings us a natural explanation of the shape of low 
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energy z/g bump observed by the MiniBooNE Collaboration. If z/2 decays into a z/i and a 
real scalar particle where z/j are neutrino mass eigenstates, it is a two-body decay where the 
energy-spectrum of the produced z/g should be discrete. It is in contrary to the observation 
where the energy spectrum of the produced z/g is continuous. Indeed, the incident z/^ beam 
has an energy distribution which can result in a natural energy spreading for the produced 
electron-neutrino, however, it demands that the shape of v^. spectrum must be similar to that 
of the incident beam. Instead, if the produced X is an unparticle, the energy spectrum of 
Vf. would naturally spread and it may be more consistent with the present measurements. The 
interactions between the unparticle and the SM particles are described in the framework of 
low energy effective theory and lead to various interesting phenomenology. There have been 
many phenomenological explorations on possible observable effects caused by unparticles 
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40[ | and much more are coming up. 
The MiniBooNE results indicate that the energy of the events of excess is about a half 
of the peak position at the energy spectrum of the muon neutrino. As discussed above, 
we suggest a decay mode V2 v\ where lA denotes the unparticle and a consequent 
transition z/^ u^+U might be observed, namely z/^ and are not physical eigenstates, but 
are that of weak interaction and can be caught by detector as an appearance of z/g at lower 
energy. As indicated in [3], the unparticle stuff with scale dimension du cannot be "seen" 
directly, it would manifest itself as a missing energy. When the scale dimension du is not 
very large, the energy spectrum of electron neutrino can fall into the allowed range of the 
MiniBooNE measurements. This process has also been considered in [2^. A transition into 
a three-body final state — z/g + z/g + Ve is also a possible process to explain the MiniBooNE 
data. The two decay modes: z/^ — z/g + W and ^ u^. + + Ve are both lepton flavor 
violating processes and can only occur via new physics mechanism beyond the SM. 

If the mechanism proposed above can explain the observed peak which depends on its 
decay width, the possible detection rate of the number of z/g is roughly 



A^o(l 



-tjr 



) X ^, 



(1) 



where iVo is the muon neutrino number, r is the lifetime which should be calculated in the 
aforementioned scenario and 17 is a detection rate which is also very small, say, 10^^° or 
even smaller, t is the flight time from the source to the detector and since the speed of 
the beam neutrino is very close to the speed of light, t ~ L/c where L is the distance and 
approximately 500 m in the MiniBooNE experiments. In addition the ratio would be further 
suppressed by the time dilation factor 7 = m/ E. Since L is only of several hundred meters, 
to make a sizable ratio which can be observed, r must not be large. We will obtain its value 
by imputing all the concerned model parameters which are fixed by fitting data of other 
experiments into our numerical computation. 

As is well-known, neutrino oscillation had been observed in the solar, atmospheric, accel- 
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erator neutrino experiments and the present theoretical studies almost completely confirm 
the MSW mechanism. The relevant mixing parameters and the mass-square differences are 
determined by fitting the data, even though the absolute values of the neutrino masses are 
still not fixed yet. While theoretically determining the parameters, possible neutrino decays 
are not taken into account seriously. How to reconcile the neutrino decay with the present 
theoretical works on the neutrino oscillation is an open question, namely, one should explore 
if there exists discrepancy between the theoretical predictions and data. Obviously if the 
decay rate is sufficiently small, one does not need to modify the present theoretical frame- 
work about the neutrino oscillation, but if the decay rate is not too small, the data fitting 
should be re-considered, therefore the MiniBooNE result indeed provides a new challenge to 



the neutrino physics and we will return to this topic in our next work 4l[ |. 



II. NEUTRINO DECAYS IN UNPARTICLE PHYSICS 

We start with a brief review of the unparticle physics. First let us consider an unparticle 
U with scale dimension du and momentum P. The unparticle momentum satisfies the 



constraint > 0. According to 18||, the unparticle stuff can be viewed as a tower of massive 
particles with mass spacing tends to zero. Scale invariance provides the most important 
constraint on the properties of unparticles. The two-point function of scalar unparticle field 
operator Ou is written as 

(0|O^(x)Oi(0)|0) = j |^^e-^^-1(0|O^(0)|P)|V(P^), (2) 
where \P) is the unparticle state with momentum P and the phase space factor is 

moump)?p{p') = Aaj{py{p^){py--\ (3) 

where 

_ IGtt^/^ T{du + l/2) 
{2T:Y^uY{du-\)V{2du)- 

For the vector unparticle field O^, we have 

(0|O^(0)|P)(P|O^(0)|0)p(p2)=A,,e(P°)e(p2)(-^^'^ + P^P7p2)(p2)'^"-^ (5) 
where the transverse condition d^0^ = is required. The Lorentz structure of unparticle 



can also be tensor [id, l32| or spinor [Ul] . In this study, we restrict our discussions to scalar 
and vector unparticles. Obviously similar analysis can be done for the tensor and spinor 
unparticles. 

About the virtual effects, the propagator of the scalar unparticle field is given as 

d^.e-(0|TO„(x)O.(0)|0) = (6) 
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and for the vector unparticle field, tlie propagator is 

The function sm{duT^) at the denominator imphes that the scale dimension du cannot be 
integers for > 1 in order to avoid singularity. The phase factor e"*^'^""^-*'^ provides a 
CP conservii ig p hase which produces peculiar interference effects in high energy scattering 
processes jo], [iB, 32] and CP violation in B decays 
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In this study, we will discuss interactions between the unparticles and neutrinos. The 
framework which describes these interactions is a low energy effective theory. For our pur- 
pose, the coupling of unparticle to neutrinos (z/^ and z/g) is given in the form as 



= T4r^/37M(1 - l^>ad^Ou + ^^/>/37m(1 - l^>aO^ + h.c. (8) 

Here, we have used the V — A type current as in the SM. The cs and cy are dimensionless 
coefficients. The z/^ and are weak eigenstates with different flavor numbers a and j3. 

As in {2^, the neutrino decay is conveniently represented in the basis of mass eigenstates 
z/j (i=l,2, we only consider two generations in this case). The interactions between unparticle 
and neutrinos are rewritten by 

^z/,7m(1 - l^^d^Ou + -^^,7m(1 - 75)^.0^ + h.c. (9) 

The relation between the coupling coefficients c's^y^ and c^'^y) can be obtained from neutrino 
mixing matrix. For a simple case, considering two neutrino mixing. 




(10) 



( cosO — sin6' 
J \ sinO cos9 

The coefficients in the mass basis are related to those in the flavor basis by 

c2?.. = cos^e (11) 



^UIU2 _ 1 l^aUll 

-SiV) ~ 2 '^S{V)- 



For a maximal mixing where 9 = 7r/4, c^^fy-, = | '^'s'('v)- '^^^ coefficients in the different basis 
differ by a constant factor. 



A. The decay of z/2 — > z^i + 

Assuming two generation neutrinos and the heavier one is 1^2 and lighter is z/i, the decay 
of z/^ — > z/e + W is realized via the transition z/2 — > z^i + W which is a typical lepton flavor 
violating process and its Feynman diagram is depicted in Fig. [H Here, it is natural to assume 
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that the basis of the interaction between unparticle and neutrino is the same as the weak 
interaction, namely v^. and z/^ are the eigenstates of the interaction. The final unparticle is 
invisible and behaves as a missing energy. The decay V2 ~^ 1^1 +U seems to be a two-body 
process. But it is different from the common case with two final particles whose momenta are 
single-valued and fixed. For the unparticle case, the energy of vi depends on the momentum 
square of unparticle which only is constrained by the condition > 0. Namely, one can 
expect that can vary within a range < < P'^^x where Pmax would be determined by 
the momenutum conservation in z/^ decay. Thus the varying P^ causes a continuous energy 
spectrum of E^^ and it is a characteristic effect of the unparticle. 

V2 




FIG. 1: The diagram for the decay oi V2 — > v\ ^lA. The double dashed lines represent the 
unparticle. 



In order to make the process realizable, the mass of should be larger than that of 
1^1, that is the so-called normal order in literature. Without losing generality, we further 
assume rriy^ ^ rriy^ and neglect rriy^ in the analysis below. The differential decay rate of 
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where the phase space factor is 



d^{p) = {2'KY8\p2-pi-q) 



2\f 2\d,j-2 



with q = P2 — Pi and the Lorentz-invariant amplitude square is 
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4 pliPi-P2), 



for the scalar unparticle and 
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for the vector unparticle. 

In the rest frame of z/2, it is straightforward to derive the differential decay rate over the 
vi energy Ei and the decay rate of V2{p2) ~^ ^lipi) 



dTs^ 
dEi 

dTv 

dEi 



^1^11^2 12 
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c 



y2i^i 12 
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du 



Sn^duid^ - 1 
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2 A ml E? 
Adu 1 
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9{m,, - 2Ei] 



(16) 
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2^1; 



(17) 



M^du{du + l){du-2) 

with du > ^ for the scalar unparticle and du > 2 for the vector unparticle. 

In this study, we concern the observation in the laboratory frame where the initial muon 
neutrino moves nearly with the speed of light. The energy of z/2 is at the order of several 
hundred MeV and is much bigger than its inertial mass, E2 ^ fn^^. The invariant amplitude 
square in the laboratory frame depends on the angle 6' between the moving directions of 
muon and electron neutrinos. Indeed, we need to do some treatments to get an analytical 
result, i.e. boost the result in the rest frame of 1^2 into the laboratory frame. The momentum 



1^2 



2E; 



-) and the momentum product 



of 1^2 is approximated by \p2\ = \lEo — = £"2(1 

P2-PI = EiE; 

0<l-cos^^'<,^^^^,. i., 
neutrino is almost parallel to that of the muon neutrino. After performing an integration 
over cos6'', we obtain the differential decay rate of 1^2(^2) — ^ ^liPi) +U{q) in the laboratory 
frame as 



+ 1 - cos^')- From > 0, the range of cos^' is determined to be 
— — ^) which means that three-momentum of the produced electron 
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for the scalar unparticle with y = ^ and du > I- For the vector unparticle, the differential 
decay rate is 



dTv 


J^1^2 


'^Adu 


f<2^ 


dEi 
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eV 



X 



1/(3 - 2y) 
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yY'^-'T{du + l) 
+ {3-7y + Ay' 
9{E2 — El). 



2Fi{l,3-du + 2-l) ^ _ ^ ^^^2^ 2^1(2, 3; du 

T{du + 2) 
2 2Fii3,3;du + A;l 



3-1) 



ridu + 3) 



(19) 



Tidu + 4) 

with du > 2 and 2Ei{a,b; c; z) is the hypergeometric function. The decay rates Ts and Fy 
can be obtained and the final results differ from Eqs. ( fT6|) and ( fTTl) by a Lorentz factor 
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B. The three-body decay of 1/2 — > z^i + i^i + J^i 

As briefly discussed in the introduction, there is another possibihty to observe a continu- 
ous energy spectrum of Ui. Now, we turn to the three-body decay of z/2 in the framework of 
unparticle: 1/2 ^ I'l + i^i + i^i- The Feynman diagram is depicted in Fig. [2] where the unpar- 
tilce serves as an intermediate agent. Because the final states have two electron neutrinos, 
there are two diagrams in the process and one needs to consider the anti-symmetrization 
of the two identical fermions. We consider only the vector unparticle part because the 
scalar unparticle contribution is proportional to the light neutrino mass and should be very 
suppressed. According to the effective interaction of neutrinos and unparticle, the decay 
amplitude of i'2{Po) i^i{Pi) + ^^1(^2) + i^iiPs) is 

M = M1 + M2, 



Ml 



^ Xl!^ ^d,e-^ u{p2h,{l - 75)^(Po)^(Pi)7m(1 - 75)^(^3) .... 

where = {du — 2)n, qi = Po — Pi and q2 = Po — P2- The amplitudes Mi and M2 represent 
contributions from Fig. [2](a) and (b), respectively. In the derivations, we have neglected 
QiQi/Qi and (3'2i?2/'?2 terms. The square of the invariant matrix element is 



\viU2Uiuii2 42 r 1 1 n2 



2^'""' 8Af/"-' 4(sind^,7r)2 



spins 



(g2)2-2dM + (g2)2-2dw 



256(po-P3)(pi-P2). (21) 




FIG. 2: The Feynman diagram for the three-body decay of z^2 ^ z^i + z^i + ^^i • 

For the three-body decays, we work in the rest frame of the muon neutrino. As will be 
shown later, the three-body decay rate of z/2 ^ z^i + z/i + ui is much smaller than that of the 
two-body decay 1/2 —>■ ui+U. This conclusion does not depend on which reference frame we 
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choose. In the rest frame of z/2, the three-body kinematics are described in terms of the final 
neutrino energies Ei, E2 by 



Po-P3 = {rriu^ - Ei - E2) , Pi ■ P2 = [Ei + E2 ^ 

qI = {Po - Pif = - 2m^^Ei, = {po - P2Y = - 2m^^E2, (22) 

Thus, the differential decay rate is 

dTy{u2 ^ z/1 + z/1 + u,) = 777^^ lllMl' dE^dE2. (23) 
where the integration range is < -Ei < and — Ei < E2 < 



Note that the simi 
have been studied in 



ar lepton flavor violating processes e +W, /i e 



15 



I22I I and their formulations are quite similar to ours. 



III. ANALYSIS ON THE CONCERNED PHENOMENOLOGY 

For the neutrino accelerator experiment, the neutrinos fly over a baseline with distance 
L before reaching at the final detectors. If neutrino decays as we suggested, the number of 
final electron neutrinos produced from muon neutrino decaying is 



N,^ = No exp 



t 



Tlab 



L m , , 



with A^o the initial muon neutrino number, riab and are neutrino life times in the laboratory 
and rest frame, respectively. In the MiniBooNE measurement, L/i? ~ 500 m/500 MeV, 
N^^ ~ 100. The ratio of neutrino life time over mass Tu/rrii, ~ -^10"^'* s/eV. When 
Nq/N,^ ~ 10^0, rjm, ~ 10-^ and when No/N,^ ~ 10^ r./m, ~ lO-^? 

At present, our knowledge on the neutrino mass is mainly obtained from the neutrino 
oscillation data. The squared mass difference of the mass eigenstates Am?- = nij — m?- 

: Am^ ~ 8 X 10"^eV^, Am^ ~ 3 x 10"^eV^. We will not use the 
LSND result (Am^ ~ leV^) since it is disfavored by other neutrino experiments and the new 
MiniBooNE measurements. From these data, we choose m2 = 50 meV as the upper limit in 
our numerical calculations. 

Firstly, we give an observation that the rate of three-body decay —>■ z/g + z/g -l- i^e is much 
smaller than two-body process —>■ Uf. + U. For an illustration, we choose the parameters 
as cs = cy = 1, A;^ = 1 TcV, du = 1.1 for scalar and du = 2.1 for vector unparticles. The 
parameters satisfy the cosmological constraints. For the vector unparticle contribution, the 
rates of two-body and three-body decays are 

ry(z/^ ^ z/e + W) = 4 X 10"^^ eV, 

Tv{Uf, ^ Ue + T^e + l^e) = lX eV. (25) 
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The three-body decay rate is more than 30 orders smaller than that of the two-body case. 
The tiny ratio is due to the very weak coupling between unparticle and neutrinos (and there 
are two such vertices for the process, see Fig. 2.) and small neutrino mass. This observation 
is analogous to the process of fi~ — e~ -|- z/^ + Ue where the decay rate is proportional to 
G'jp and m^. If we only use the three-body decay, a life time of neutrino which is so long 
that muon neutrino will never decay and the bump observed in the experiment cannot be 
explained by the neutrino decays at all. Thus, we can safely neglect the contributions from 
the three-body decays and approximate Ti, = r(i/^ — z/g -|- W). 

Secondly, we discuss the constraints of unparticle parameters from neutrino decays. As 
discussed above, the MiniBooNE experiments put a bound for neutrino life time in the rest 
frame T^/m,y ~ 10~^ s/eV. We take this bound for our analysis and discuss three possibilities. 
In order for the illustration, we are restricted in the case of vector unparticle. (1) We fix 
du = 2.1, Au = I TeV, and constrain cy by cy > 10^^. The coupling constants are found 
to be much larger than the order 1. (2) We fix C5 = cy = 1, Ak = 1 TeV, and constrain 
< du — 2 < 10^"^. The scale dimension will be nearly equal to 2. (3) We fix cs = cy = 1, 
du = 2.1, and constrain Ay < 10^^ TeV which is obviously impossible. Thus, if neutrino 
decays as suggested, the unparticle parameters have to fall into a very unnatural space. On 
the opposite side, if the unparticle parameters are chosen in a reasonable space, the neutrino 
life time will be so long that they will not decay when flying over the distance L ~ 500m in 
the MiniBooNE. 

Thirdly, we discuss the relative energy spectrum of final electron neutrino. Fig. [3] plots 
the initial z/^ energy spectrum. The distribution is a quasi-Gaussian function, which peaks 
around 700 MeV. If the muon neutrino decays to electron neutrino and a conventional 
particle, the energy spectrum of electron neutrino has the same distribution as that the 
muon neutrino beam. From the data of excess events plotted Fig. HI obviously, the z/g 
energy spectrum is not consistent with the data, no matter for low energy or high energy. 
The experimental data show that the excess of events is a decreasing function rather 
than a Gaussian distribution in the energy range 0.3 < E^^^ < 1.0 GeV. This excludes the 
neutrino decays with conventional particles. Unparticle is not a regular particle and has 
no a fixed mass. The energy distribution of z/g with unparticle being in the final state is 
different from the initial z/^ spectrum and the final energy distribution of z/g depends on 
both effects. Combining the z/^ energy spectrum (energy spreading of the muon neutrino 
beam) and the differential ratios of neutrino decay given in Eqs. fll8|19l) . the final electron 
neutrino energy distribution is depicted in Figs. H] and [5] for scalar and vector unparticles, 
respectively. Within the non-oscillation region, i.e., at low energy 0.3 < E^^ < 0.45 GeV, 
the theory prediction is well consistent with the experimental data. For the energy range 
0.5 < Ei,^ < 1.05 GeV, the theory prediction is slightly larger than the data. It is noted 
that in the above figures, only the relative size is estimated, the absolute magnitude is very 
small. 
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FIG. 3: The energy spectrum for the beam in the laboratory frame. 
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FIG. 4: The energy spectrum for the decay of — >^ t'e +W with scalar unparticle in the laboratory 
frame where du = 1-1, A^^ = 1 TeV and cs = cy = 1- 

IV. DISCUSSIONS AND CONCLUSIONS 

Motivated by the new measurement of the MiniBooNE Collaboration, which observed an 
excess of electron-like events at low energy, it motivates us to search for a possible mechanism 
beyond the SM to explain this phenomenon, so the first idea which hits our mind is that U2 
might decay into z/i accompanied by some other very light products. We should testify if 
this scenario can produce results which are theoretically plausible and can explain the data. 

There may be several possible modes, the first one is 1/2 — > t^i + z/i + z/i which is a three- 
body decay, the second one is 1^2 — > i^i + a where a is a single boson particle, for example 
an axion etc., and the third one is 1/2 ^ 1^1+ U where U represents an unparticle. All the 
three possibilities cannot be realized in the framework of the SM, so new physics beyond 
the SM is necessary. The first one was numerically estimated in this work and our results 
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FIG. 5: The energy spectrum for the decay of — > +IA with vector unparticle where du = 2.05, 
A.U = 1 TeV and cs = cy = 1. 

indicate that the decay rate determined by the three-body decay mode is too small and 
is ruled out immediately. The second mode is a two-body decay, therefore the spectrum 
of the electron neutrinos is discrete and it is not consistent with the measurement of the 
MiniBooNE. Even though we consider the energy spreading of the incident muon neutrino 
beam, the shape of the resultant electron neutrino bump cannot be well understood in this 
scenario. Therefore the third candidate is the most favorable. In this work, we work out the 
formulations of neutrino decays within the framework of unparticle physics. The formulations 
in the laboratory frame are given for the first time. 

The smallness of the decay width given by our numerical results indicates that the un- 
particle scenario may not explain the excess of electron neutrinos at low energy. The life 
time predicted in the unparticle model is qualitatively consistent with the the cosmological 
constraint 43|] which is about 10^^ sec. By eq. (1), we know the suppression of exp(— 
with t ~ 10~^ sec, would kill any possibility of observing a decay event. The reasons are: 
(1) very tiny neutrino mass (2) very weak interactions between the unparticle and neutrino. 
Since our numerical results are consistent with the cosmology constraints and the results by 



other authors 26|, we can be convinced that the calculation is right, but the proposal does 
not work here. 

Thus in this work, we definitely obtain a negative conclusion that the peak of electron- 
neutrino at lower energy observed by the MiniMoone collaboration cannot be explained 
by the neutrino decay. On other aspect, the phenomenon is there and demands theoretical 
explanations, so that we propose another scenario which might overcome the aforementioned 
restrictions which forbid the appearance of electron neutrinos to appear at low energy for 
the MiniBooNE experiments. We will present the scenario in our next work. 
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